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Abstract. Let p(ui) denote the probability that four random circular 
caps of angular radius 70° < oj < 90° cover the unit sphere S 2 . An exact 
expression for p(uj) is unknown. We give nontrivial lower bounds for p(u>) 
when u> > 84°; no improvement on the inequality p(uj) > 0 for ui < 84° is yet 
feasible. A dual problem involving randomly inscribed well-centered tetrahedra 
is also examined. 


Let S 2 denote the two-dimensional sphere in M 3 of unit radius. Let Xi, X 2 , A 3 , 
X 4 be four points that are independent and uniformly distributed over the surface of 
S 2 . Each X) is taken to be the center of a circular cap of angular radius uj\ 

{Y G S 2 : arccos (Xj ■ Y) < cu } , 0 < lo < it 


where A" • Y is the usual inner product between vectors in M 3 . A point, hemisphere 
and S 2 itself are the outcomes if u — 0, 7t/2 and n. The sphere is said to be covered if 
each point of S 2 belongs to at least one cap. It is known that the coverage probability 
p(u)) satisfies [L3| |3j 


where 


p(u) 


f 0 if 0 < oj < oj 0 , 

\ 1 — 2 cos 6 (|) [8 — 9 cos 2 (|)] if f < u < vr 

co 0 = arccos (|) = 1.23095941... » 70.53°. 


Finding an exact expression for p(oj) in the interval u 0 < uo < n/2 remains an open 
problem. Early works in this area include [Ml E]; computer simulations appear in 
PSl E2] and historical overviews in da ei is]. All such papers deal with n circular 
caps, where n is often large. Our restriction to the case n — 4 is somewhat unusual, 
as far as is known. 

Gilbert [5] gave the following bounds on p[pj): 


2 

0 < p{pj) < ^(fc)(-l) fc (l - k sin 2 (f )) 4 = 1 - 2 cos 8 (|) + cos 4 (tu) 

k =0 
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(his left-hand side is negative for n = 4 and thus replaced by zero). The upper bound 
is weak but not relevant to us. We will focus on the lower bound. Our improvement, 
based on elaborate formulation, occurs only when oj > 84.25°. For example, if 
oj = 88°, then 0.0765 < p{oj) < 0.8567 via our work and p{oj) ~ 0.08 via computer 
simulation. 

Dual to the coverage problem is the following. Given A, B , C , D on S 2 , there 
is a unique circular cap of minimal angular radius 9 m ; n containing all four points. If 
A, B, C , D are independent and uniform on S 2 , then the cumulative distribution 
function $ for the random variable 9 min satisfies PHIS] 

1 — <E>(7T — oj) — p{pj), o < OJ < IT. 


Define a distribution 

HO = 8 {m - |) = 8P{f < 9 min < e} , f < £ < VT. 

Calculating F(£) is hopeless; therefore we study a different variable 9 ahc with distribu¬ 
tion G that is computationally feasible and satisfies cG(£) > F(£) for some constant 
c > 0. 

Any three points from {A, B, C, D} almost surely form both a chordal triangle 
(with sides as straight lines through the interior of S 2 ) and a spherical triangle (with 
sides as great circles on the surface of S 2 ). 

The chordal triangle ABC determines a plane that partitions S 2 into two com¬ 
plementary caps. Let ABCD denote the tetrahedron inscribed in S 2 with base ABC 
and apex D. Such a tetrahedron is called well-centered if it contains the origin 

m- 

The hypothesis 7t/2 < 9 min implies that ABCD is well-centered and, without loss 
of generality, ABC is acute. Further, D is in the larger of the two complementary 
caps and the angular radius 9 &bc of the cap satisfies tt/2 < # abc . The converse is not 
true. Our main objective is to find the conditional probability density function g for 
9 a be and the associated distribution G. Since ABC is one of four possible choices of 
base, we have a crude estimate 4G(£) > F(£). 

While finding g , we compute the probability 

k = P {ABCD is well-centered & ABC is acute} = 0.10191818... 

which extends a well-known result P {ABCD is well-centered} = 1/8 from [2T| [20 , 
mmm and another result P {ABC is acute} = 1/2 from [13]. Details underlying 
the acuteness probability are given in the next section. The expected number of 
acute faces of ABCD, given well-centeredness, is 32k; a refined estimate is hence 
(3.26138193.,.)G(0 > F(£). 
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1. Random Triangles of Circumradius r 
For comparison’s sake, consider triangles formed by selecting three independent points 
on the circle r S 1 , where the radius r is fixed. The bivariate density for two arbitrary 
angles a, f3 is [3] 

f 2 

I — if0<a<7r, 0</3<7t and a + (3 < n, 

■\ n 2 

I 0 otherwise 


thus the acuteness probability for such triangles is 1/4. The bivariate density for 
two arbitrary sides a, b is 


4 1 


0 


if 0 < — < 2 and 0 < - < 2, 

r r 

otherwise; 


the sides a, b are independent despite the fact that angles a, (3 are dependent and 
a = 2rsin(a), b = 2rsin(/3). 

Consider now chordal triangles formed by selecting three independent points on 
the sphere S 2 . Given triangle ABC , let r denote the radius of the unique circle passing 
through vertices A, B, C. Unlike before, r is a random variable. Intuition might 
suggest that the density of a, f3 and of a/r, b/r should be similar to the preceding. 
In fact, however, the bivariate density of a, [3 is 


— sin(a') sin(/3) sin(a + f3) 
37T 
0 


if0<a<7r, 0 < /3 < 7T and a + (3 < tt, 
otherwise 


thus the acuteness probability for such triangles is 1/2. The trivariate density of a, 
b, r is 



for 0 < a < 2r and 0 < b < 2r. This is a (2/3, l/3)-weighted mixture of densities. 
The portion coming from acute chordal triangles is 


ab iU-MT + R-M;) 2 

6rr4 


a b 
Qirr A 


a 

y/4 r 2 — a 2 


+ 


b 


\J4 r 2 — b 2 

























Covering a Sphere with Four Random Circular Caps 


4 


for 0 < a < 2r, 0 < b < 2r and a 2 + b 2 > 4r 2 . We call this portion 5(a,b,9) for 
future reference, where r = sin($). 

Proof of the above follows via a Jacobian determinant calculation based on formula 
(4.18) of [T2] : see also Theorem 3.2 of 0 - This is one of several crucial contributions 
by Miles to our study. 

2. Conditional Probabilities and Angular Radius 
Let £ denote the event {ABCD is well-centered & ABC is acute}. Following the 
proof of Theorem 3.3 of |13j . we have 

P{0 < #abc < 0 + d6} ■ F{£ | 9} = P{0 < 0 abc <9 + d9\£}- P{£} 
where t/2 < 6 < n. Formula (3.3) of [13J yields 

P{ 6 > < # abc < 9 + d9} = | sin 3 ($) d9 

and finding 

P {0 <1 0 a bc ^ 9 + d9\£} = g(9) d9 

is our main objective. We turn attention to P{£ | 9}. 

Fix a chordal triangle ABC for consideration. Let A denote the spherical triangle 
with vertices A, B, C. The set of all points D for which the tetrahedron ABCD is 
well-centered is A', the spherical triangle with vertices A', B', C antipodal to A, B, 
C. This pictured in Figure 2 of [TH] and is stated in [I]. Note that the area of A' is 
equal to the area of A. 

Letting ABC vary over all acute chordal triangles of fixed angular radius 7 r /2 < 
9 < 7 r, the random quantity area(A) can be as large as —6 arccot (\/3cos0) —i r. This 
is the area of an equilateral spherical triangle with vertices on a circle of radius sin( 0 ) 
[ 2 ], A formula for the mean area, as opposed to the maximum area, turns out to be 
quite complicated. We derive this in the next section. 

3. Kolmogorov-Robbins Theorem 

The desired probability P{£ | 9} is the expected value of area(A)/( 47 r), where random 
spherical triangles A correspond to acute chordal triangles ABC of fixed angular 
radius 7t/2 < 9 < tt [ 10 } [16} 8 j ITTj . 

Let r = sin(0) for convenience. If angles a, /3 of an acute chordal triangle are 
known, then the remaining angle 7 satisfies 

cos( 7 ) = cos (t — arcsin (^) — arcsin (^)) 
ab — \/4r 2 — a 2 \/4r 2 — b 2 
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The corresponding spherical triangle has angle 7 given by [T5j 


7 = arccos 

= arccos 


4 cos( 7 ) — a b 
( 1 - 


- a 2 \J 4 — b 2 
r 2 )a b — \/4r 2 — a 2 \/4r 2 — b 2 
r 2 \/4 — a 2 \/4 — b 2 


Call this expression A (a,b,9). While area(A) = a + /3 + 7 — tt and analogous 
expressions for a, /3 are possible, it is simpler to employ E(area(A)) = 3E(y) — tt. 

This yields 

2 r 2 r 

P{£ | 9} = I I — (3\(a,b,9) — tt) 6{a,b,9) dadb 

0 yj 4^*2 - Jj‘2 

and 

7r 

P{£} = J ^ sin 3 (6 ) ) P{£ | 9} d9 = 0.10191818... = k: 

7t/2 

therefore we conclude that 
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G(9) = 



7r/2 


[ — sin 3 (t) P{£ I t} dt. 
J 2k 


In Figure 1, a plot of the density g(9) is superimposed on a histogram of simulated 
6 >abc-values, 9 e [7 t/2,7t], based on 10 6 random well-centered tetrahedra ABCD with 
acute faces ABC. The fit is excellent. First and second moments are also indicated. 
All required integrals were evaluated via numerical methods because no simplification 
via computer algebra seems likely. 

Of course, our real interest is in 9 min , not 9 a b c . Define a function 


m 


24sin 5 (f) cos (f) [2 - 3 sin 2 (f)] if 0 < 9 < f, 
|(32 n)g(9) if | < 9 < n. 


In Figure 2, a plot of 9 ) is superimposed on a histogram of simulated 0 min -values, 
9 G [ 0 , 7 r], based on 10 6 random (arbitrary) tetrahedra ABCD. The curve ifr is 
represented by the solid line and the fit is excellent on [0,7t/2], The area under tJt 
over [0,7t/2] is 7/8; the area under ip over [7t/2,7t] is 4k = 0.40767274.... It visibly 
dominates the (unknown) density function <p for 9 min . 











Density 
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Experimental mean is 1.81225; mean square is 3.32669 
Theoretical mean is 1.81212; mean square is 3.32619 



1.75 2.00 2.25 2.50 2.75 3.00 


Angular Radius theta.abc 


Figure 1: Experimental histogram of # a b c -values and theoretical fit. 



Density 
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Experimental mean is 1.28948; mean square is 1.72901 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Angular Radius theta.min 


Figure 2: Experimental histogram of 0 m i n -values and theoretical fit/upper bounds. 
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4. Bounds on Coverage Probability 

Define a function 


q(oj) = | (1 — 32 k G(tt — uj)) , 0 < u < | 

then q(oj) is a nontrivial lower bound for p(cv) when q(uj) > 0, that is, when uj > 
84.25°. As discussed earlier, q( 88°) = 0.0765. 

Our method fails to improve upon the inequality p(u>) > 0 for u> < 84.25°. We 
assess the effect of more closely approximating the density ip. Let 


e 0 = 7T-uj 0 = arccos (-§) = 1.91063323... w 109.47°. 


Assuming the function 1 + p is logarithmically convex on [7 t/2,0o]> it follows that 


V’lcvW 



e 0 -e 

Bq—tv/2 


- 1 


dominates ip. This is suggested in Figure 2, where the curve i/j \ cv is represented by 
the dashed line. Integrating ifi\ cv from 7 t/ 2 to 9 , we obtain \Eh cv , which in turn is used 
to define 

gicvM = I - Ti cv (7T - a;), 0 < u < f. 

This is a tighter bound than the preceding; for example, gi cv (88°) = 0.0766. It is 
nontrivial when q\ cv (co) > 0, that is, when lu > 83.90°. Thus the effect is only slight. 
A rigorous proof that log convexity holds, however, is not known. 

A sharper upper bound for p(u>) than that provided in [5] is also desired. 


5. Addendum 

The preceding text was written in 2011. A question, “Does k possess a closed-form 
expression?” appeared in [23]; an affirmative answer 

11 1 
K ~ 96 _ 8 ^ 

appeared in [24] with a different line of reasoning. We acknowledge that the preceding 
text is heuristic and non-rigorous. Sample open issues include: Must a well-centered 
tetrahedron ABCD possess at least one acute-angled triangular face? (Simulation 
suggests, in fact, that there are at least two.) If all faces ABC , ABD, ACD, BCD 
are acute, then 

dnin mill {ffibc? ^abd; Cicd • 9\ iC( \ | 

and the inequality 4G(£) > F(£) becomes clear. It is possible, however, that only 
two or three faces are acute. Let N e {2,3,4} be the number of acute faces, given 
a random well-centered tetrahedron ABCD. What precisely is the distribution of 
N? All we know is E(N) = 32k. How is the inequality (32 k)G(0 > F(0 proved? 
(The apparent continuity of ip {9) at 9 = 7t/ 2 gives only an impression that this may 
be correct.) Discussion would be appreciated. 
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